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RESTRAINED AGAINST ROTATION 
By S. B. Batdorf and John C. Boubolt 


SUMMARY 


An exact solution and a closely concurring approxi- 
mate energy solution are given for the buckling of an 
infinitely long flat plate under combined shear and 
transverse direct stress with edges elastically restrained 
against rotation. It was found that an appreciable frac- 
tion of the critical stress in pure shear may be applied 
to the plate without any reduction in the transverse 
compressive stress necessary to produce buckling. An 
interaction formula in general use was shown to be 
decidedly conservative for the range in which it is .sup- 
posed to apply. 


INTRODUCTION 


In the design of- stressed-skin structures, considera- 
tion must sometimes be given to the critical stresses for 
a sheet under a combination of shear and direct stress. 

The upper surface of a wing in normal flight, for example, 
is subjected to combined shear and compressive stress, 
and the lower surface is subjected to combined shear and 
tensile stress. The upper surface may then buckle at a 
lower compressive or shear stress than if either stress 
were acting alone. The critical shear stress for the 
lower surface may be increased by the presence of the 
tensile stress. 
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If the wing has closely spaced chordwise stiffeners 
the skin between two adjacent stiffeners may be regarded 
as a long sheet slightly curved in the longitudinal or 
chordwise direction and straight in the transverse or 
spanwise direction loaded in shear and- transverse direct 
stress. A conservative preliminary estimate of the 
critical stresses may be obtained if the sheet is con- 
sidered to be flat and infinitely long. In the present 
paper, the critical stresses are computed for an infin- 
itely long flat plate loaded as indicated in figure 1(a) 
The corresponding idealization of the case of a wing 
with spanwise stiffeners, shown in figure 1(b), was 
treated in reference 1. 


Ls Jbr — 1 — 1 sasJLfei t •-*- 1 =& J l 1 JLsJ 


1 ! 




F^T T^T^' I ■ t”'*’ f i ^ f p^T f • f 

(a) Type of loading problem solved in present paper, 



(b) Type of loading problem solved in reference 1. 


Figure l.-r Buckling of an Infinitely long plate under 

combined loads. 
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CONVENTIONAL INTERACTION FORMULAS 


^For buckling of structures under combined loading 
conditions, no general theory has been developed that 
is applicable to all cases. Stress ratios (reference 2), 
however, provide a convenient method of representing 
such conditions. For example, the ratio of the shear 
stress actually present in a structure to the critical 
shear stress of the structure when no other stresses are 
present may be called the shear-stress ratio. Stress 
ratios may similarly be defined for each type of stress 
occurring in the structure. 

It is generally assumed that equations of the type 


R]_P 


Rs q 


+ R-! 


= 1 


( 1 ) 


may be used to express the buckling conditions in the 
case of combined loading (reference 3, p. 1 - 18) . In 
equation (1), Rj , Rg, and R^ are stress ratios and 

p, q, and r are exponents chosen to fit the known 
results. (All symbols are defined in appendix A.) 

Such a formula gives the correct results" when only one 
type cf' loading is present and has the further advantage 
of being nondimensional . Equation (1) implies, more- 
over, that the presence of any positive fraction of the 
critical stress of one type reduces the amount of 
another type of stress required to produce buckles; this 
implication appears reasonable and has been proved true 
in some cases (references 1, 2, and 4) . 

In reference 2 the following interaction formula is 
given for an infinitely long plate with clamped edges 
loaded in shear and longitudinal compression: 


Rg 1 * 5 + R c = 1 


( 2 ) 


where R g is shear-stress ratio and R c is longitudinal 

direct- stress ratio. The same formula is recommended in 
reference 3 for general use for the buckling of any flat 
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rectangular plate, regardless of the direction of com- 
pression and the degree of edge restraint. 

Later theoretical work (reference 1) shows that, to 
a high degree of accuracy, for an infinitely long plate 
with any degree of edge restraint loaded in shear and 
longitudinal compression 


R s 2 + R c ■= 1 (3) 


The same formula was found in reference 5 to be appli- 
cable to simply supported rectangular plates of aspect 
ratios 0.5, 1, and 2 ; and the conclusion was drawn that 
interaction curves in stress-ratio form are practically 
independent of the dimensions of the plate. 

The present analysis, however, indicates that the 
buckling of an infinitely long plate loaded in shear 
and transverse compression (fig. 1(a)) is not adequately 
represented either by equation (2) or (3) or by any 
formula of the type of equation (1) . Two independent 
theoretical solutions to this buckling problem are given 
in appendixes. B and C. Appendix B contains the exact 
solution of the differential equation of equilibrium, 
and appendix C contains an energy solution leading 
directly to an interaction formula. This energy solu- 
tion, which gives . approximate values only, was made to 
obtain an initial quick survey of the problem and to 
provide a check on the results of the exact solution. 
Approximate interaction formulas in substantial agreement 
with these results were given for the cases of simply 
supported and clamped edges in reference 6. 


RESULTS AND DISCUSSION 


In figure 2, curves are given that indicate the 
critical combinations of shear and. transverse direct 
stress for an infinitely long plate with edges elasti- 
cally restrained against rotation. These curves are 
computed from the exact solution presented in appendix B. 
The d.egree of edge restraint is denoted by e, which 
is defined in appendix B in such a way that zero edge 
restraint corresponds to simply supported edges and 



NACA ARR No. L4L14 


5 


infinite edge restraint indicates clamped edges. A 
similar set of curves is given in terms of stress ratios 
in figure 3. The numerical values used to plot figures 2 
and 3, together with the values found by the energy 
solution, are given in table I. 

The most striking feature of these results is that 
an appreciable fraction of the critical stress in pure 
shear can evidently be applied to the. plate without any 
reduction in the compressive stress necessary to produce 
buckling. (See fig. 3.) This fraction varies from 
about one- third to more than one-half, depending on the 
degree of restraint. At shear stresses higher than those 
corresponding to this fraction, the compressive stress 
required to produce buckling is reduced by the presence 
of shear. The result that the compressive buckling 
stress is entirely unaffected by the presence of a con- 
siderable amount of shear is probably peculiar to infin- 
itely long, plates. It is to be expected, however, that 
this result will be closely approached in the case of 
long finite plates. 

In figure 4 a comparison is made between the exact 
solutions and the interaction formulas of equations (2) 
and (3). Equation (2), which is the. interaction formula 
in general use, is seen to be decidedly conservative. 


CONCLUSIONS 


The exact solution of the differential equation for 
the buckling of an infinitely long flat plate under com- 
bined shear and transverse direct stress with edges 
elastically restrained against rotation indicates the 
following: 

1. An infinitely long flat plate may be loaded with 
an appreciable fraction of 'its critical stress in pure 
shear without causing any reduction in the transverse 
compressive stress necessary to produce buckling. 

2. An interaction formula in general use for rec- 
tangular plates in .combined shear and compression is 
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decidedly conservative when applied to an infinitely long 
plate in shear and transverse compression. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX A 


SYMBOLS 


C 1» C 2 
• D 


E 

N 

N 


y 

xy 


S 

T c 

T s 

V 1 

v 2 

Y 


c 

k c , k g 


functions of edge restraint coefficient 
given in appendix C 

flexural stiffness of plate per unit length, 

Et® | 

in- lb . o\ ; 

[ 12 (l - K^J 

elastic modulus of material, P s - 

compressive force per unit length, lb/m. 

shearing force per unit length of plate, 
lb/ in. 

_ . . ,.n p o q -nef inch of restraining 

rotational stiifness P®* ; u \ 0 / radian 

member at edge ox P late » 1J/ 

work done by compressive force per half wave 
length, in- lb 

work done by shear force per half wave length, 
in-lb 

• „ r-iTfri-p ner half wave length, 
strain energy m plare per 

in-lb 

P j,,p restraint, per half wave 
strain energy in ea& e resudil ’ - 

length, in- lb 

function of y associated with deflection of 
plate during buckling 

width of plate, in. 

width of plate in oblique coordinate system of 
reference 1, m. 

function of a, 3, and X 

critical compressive and shear-stress coeffl- 
cients, respectively 
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m root of a characteristic equation of appendix B 

fl,f 2 >f .3 functions of restraint coefficient e given, 
in appendix G 

t thickness of plate, in. 

irb 


w 


X 

y 

a,p 

e 

Y 

\ 


V- 


a 



T 


0 


displacement of buckled plate from original 
position 

amplitude of assumed wave form of buckle 
longitudinal coordinate of plate 

t r ans ve r se_.coor dina te_. of pi. ate 

functions of X, y, and k c 

nondimensional coefficient of edge restraint 

one of two parameters determining buckle form 

one of two parameters determining buckle form 
(half wave length of buckle, in.) 

Poisson's ratio 

direct stress 

transverse direct stress, psi 
shear stress, psi 

angle between buckle node and y-axis 


6 = tan 0 

^1 , ^2»^3 stress ratios 
R g shear-stress ratio 

R c longitudinal direct-stress ratio 

cr critical (used as subscript) 
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APPENDIX E 


SOLUTION BY DIFFERENTIAL EQUATION 


Statement of problem.- The exact solution for the 
critical stress at which buckling occurs in a flat 
rectangular plate subjected, to combined shear and com- 
pression in its cwn plane may be obtained by solving 
the differential equation that expresses the equilibrium 
of the buckled plate. The plate is assumed to be infin- 
itely long, and equal elastic restraints against rotation 
are assumed to be present along the two edges of the 
plate . 


Differential equation.- Figure 5 shows the coordinate 
system used. The differential equation for equilibrium 
of a flat plate under shear and transverse direct stress 
is (from reference 7) 


v ,* 4 


2 + 
\ 2s 2 

6x cy- 


6 4 w\ + 

2tt d2 * 

+ V ^ = 

0 

(Bl) 

e>y V 

ox oy 

J \ 2 



write 

Oy and 

t in terms 

of 

the 


dimensionless buckling coefficients k c and k p oy 
means of the relations 


a 


v 


k n tr c D 

1st 


k 0 TT 2 D 


(B2) 


b 2 t 


Substitution of the expressions for and r 

from equations (B2) in equation (Bl) gives 
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c 4 w + g d 4 w + b 4 w + " c g d^w + d 2 w _ 

dx 4 dx 2 dy^ dy 4 b 2 bxby ■ ^ 2 ^y2 


Solution of differential equation .- If the plate is 
infinitely long in the x-direction, all displacements 
must be periodic in x and the deflection surface may 
be taken in the form 


w = Ye 


. TTX 

1 ir 


( B4) 


where Y is a function of y only and \ is the half 
wave length- of - the— buckhe-s in - the x-drrec tlrorr. — 

Substitution of the expression for w from equa- 
tion ( B4) in the differential equation (B3) gives the 
following as the equation that determines Y: 


a 4 y 

ay 4 


+ ,'Yh . <S + o (B5) 


X 


2 ) j„2 


dy 


b 2 \ d y v 


A solution of equation (B5) is 


img 
Y = e b 


where m is a root of the characteristic equation 


m 4 + 


#) S - 


TT 2 k, 


+ (^) 4 = o 


( B6) 


Except for the substitution of 

J 2b f 


2 0 


Zk f _ tt 2 

ru ) 


TT^k r 


for 2j^- ) , equation (B6) is identical with equa- 
tion (A-6) of appendix A of reference 8, in which 
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equation (53) of this appendix was solved with k c = 0. 

“With this change, all the results obtained in appendix A 
of 'reference 8 are applicable here. The stability 
criterion for combined compression and shear is therefore 
the same in form as that for shear alone, given by equa- 
tion (A-19) of reference 8, which is 


2aS (4y 2 - -^-^(cosh 2a cos 2(3 - cos 4y) - [V( p 2 - a 2 ) 


-(p 2 ♦ a 2 f - (4Y 2 - p 2 + a 2 ) 


sinh 2a sin 2p 


+ eja(4y 2 + a 2 + p 2 ) cosh 2a sin 2p 

+ p (4y 2 - a 2 - p 2 ) sinh 2a cos 2p - 4apy sin 4yJ = 0 (B7) 


The relation between k s and a., p, ana y is 
also the same in form as that in equation (A-23) of . 
reference 8, namely. 


_ 8y(q^ + p 2 ) 


IT 


TTb 

X 


(B8) 


In the present report, however, a and p have the 
following values: 



c 


+ 


l/TTb\ 2 

4 \ X J 



8 


k 


c 


vfnere 
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As in reference 8, the restraint coefficient t is 
defined herein by the relation 



where S is the ratio of a sinusoidally applied moment 
to the resulting sinusoidally distributed rotation of 
the restraining element measured in radians. 

Evaluation of k R corresponding to a selected 
value of k c . - The procedure for evaluating k s , after 
values of k c and e have been chosen, is as follows: 

A value of bA is selected; a series of values of y 
are assumed until one is found that, together with the 
corresponding values of and 6 computed from equa- 

- tion (B9) , -sa-tisf ien-equabi-on- (-B7)-; — k- s - -i-s -then— c-om- — 

putsd from equation (R6) . Another value of bA is 
selected; a new set of values of y, a, and p is 
found that satisfies the stability criterion; and a new 
value of kg is computed. The entire process is 
repeated until the minimum value of kg can be found 
from a plot of k c against b/X, When £ is a func- 
tion of b/X, e muse be reevaluated each time a 
different value of b/X is selected. The minimum value 
of k s and the chosen value of kc, when inserted in 
equations (B2) , give a critical combination of shear 
and direct stress. 

Evaluation of k s when k c has value corresponding 
to buclEI lir^^’Tullir 'c'oT umn r^" One critical combinati o n 
of shear ana compression is simply k s = 0 and k c 

equals the value corresponding to buckling as an Euler 
column. The curves giving critical stress combinations, 
however, did not appear to be approaching this point as 
their construction progressed. It was therefore neces- 
sary to determine whether values of k s other than 
zero are critical when the Euler compressive stress is 
reached. The determination of k s when k c reaches 

the value at which the plate buckles as an Euler column 
requires special treatment, because kg, given by 
equation (38), becomes indeterminate when the wave 
length becomes infinite ag suggested by the energy 
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solution. The result that X becomes infinite when k c 
takes its Euler value is readily checked from equa- 
tion (Cll) for the special case of e = 0; for this 

case p = q = r = and (^c) cr = Prom equation (B8) 

it is clear that, if k s is to remain finite when the 
wave length approaches infinity, either 

Y — > 0 ( a) 

or 

a 2 + p 2 — >0 (b) 


For case (a), when e = 0, it follows from equa- 
tions (B9) that, to small quantities of the second order. 


where 



- “(3 y 2 + 2u 2 )] 


tr 

[3 = Y 


> 


(BIO) 


u = 


rrb 

2 \ 


If the values of a and p from equations (310) are 
substituted in equation (B7) and the resulting equation 
is expanded, with only the lowest powers of u and y 
retained., 


p 4tt 2 u 2 

Y- - : 

64 - 6tf' 


(Bll) 


Substitution of the values of a and p. from equa- 
tions (BIO) and (Bll) in equation (B8) gives as the 
final result for e = 0 


a 


k« 


= + 


2tt 


V64 - 6tr 2 
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Case (b) can be analyzed by a similar method, but 
the analysis is quite complicated because terms of third 
order must be retained. For e = 0 and e = «, case (a) 
and case (b) were found to lead to exactly the same 
result for k s . A value of k s ' other than zero when 

k c takes its Euler value may be found in the same manner 

for other values of edge restraint. For any value of the 
restraint coefficient e. 


f €TT \/ k c + 2rr \/ k c + — ^Lr^sin rr\/k c 

_V W k c/ 

(v - - 2€ - v) cos ■ -~-( c ° s n \i~c - 1) 

' “ / n'-k- 


[v~ T5 V k c - 6t W k c - 2etT\/kZ + \J k c + ^~ \/ k c^ sin Tt\f~k c 


(B12) 


, ( 3tf 2 , , etr^kc 3 2 . - 

VT ° ~T“r )cos rrvjk, 


+ + 8e + 2e^ + — - ^(cos rr \j k c - l) 


2rrk„ 


where k c has the value corresponding to Euler column 

buckling, at this value of e. • The relationship between 
this Euler value of kc and the corresponding e is - 
given by the equation (from reference 9) 
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V.hen k c reaches the Euher value , the critical shear- 
stress coefficient k s can therefore be either 0 or 

the value given by equation (B12) . The conclusion that 
k s can also have any value between these limits is 

plausible on the basis of the following physical con- 
siderations. The shear stress does no work during 
buckling when the stress condition is such that the 
plate buckles with an infinite wave length. The effect 
of shear, furthermore, is to reduce the wave length to 
a value of the order of the width of the plate. The 
wave length at the time buckling occurs, is infinite, 
however, when the plate is either in pure compression or 
at the value of k s satisfying equation (3312) . This 

fact means that, for values, of k s between 0 and that 

given by equation (B12) . the shear stress is not great 
enough to force buckling in short waves and therefore 
does not assist in producing buckles. In this range of 
shear stress the compressive stress necessary to produce 
buckling is, consequently, the Euler stress. 
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APPENDIX C 


SOLUTION BY ENERGY METHOD FOR EDGE RESTRAINT 
INDEPENDENT OF WAVE LEN GTH 


The critical stress is determined on the basis of 
the principle that the elastic-strain energy stored in a 
structure during buckling is equal to the work done by 
the applied loads during buckling. If the structure 
under consideration is an infinitely long plate under 
combined shear and edge compression with edges elasti- 
cally restrained against rotation, this equality may 
be written 


m 

A c- 


+ - T s-=- V l- + -V 


- -(-C.1-) - 


■ In reference 1 an energy solution was given for 
the type of loading shown in figure 1(b) . The deflection 
function used in reference 1 is also suitable for appli- 
cation to the solution of the type of loading shown in 
figure 1(a), which Is the loading considered in the 
present paper. The values for T g , Vp, and Vo may 

accordingly be taken directly from reference 1, but T„. 

must be recomputed to apply to the case of transverse 
compressive stress. 


The following substitutions are used to transform 
the energy expressions from the oblique coordinates of 
reference 1 to the rectangular coordinates used in the 
present paper (fig. 6): 


For brevity 


Reference 1 
y/'Oq = 

bp cos cp = 
the following notat 
tan cp = 


Present paper 

y/b ) 

b / 

ion is also adopted 


(C2) 
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By use of equations (C2) , the expressions from 
reference 1 that are used in the present paper may be 
rewritten as follows: 



and € is the restraint coefficient defined in appendix B. 

The work done by the compressive force per half wave 
length may be written 



As in reference 1, the assumed deflection function is 
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cos 


TT 

\ 


(x + 0y) 


(C7) 


When the expression for w from equation (C7) is 

substituted in (C6) and the indicated operations are 
performed, 


T = 
A c 


w, 


2 


N v rrb 

k? 

4\ * 



2 + 



( C 8) 


Values from equations (03), (04), (05), and (08) 
are now substituted in the buckling equation (Cl) . The 
use of the equations 


Nxy' 


_ Q — 

TT D 


= k 


s b 2 


N.,, 


TT^D 


c 


eliminates N X y and Ny. The resulting equation gives 

the critical combination of stresses and may be written 
as 



This equation shows that, for a selected value of k c , 

the critical shear stress depends upon the wave length 
and the angle of the buckle. Since a structure buckles 
at the lowest stress at which instability can occur, 
k g is minimized. with respect to wave length and angle 
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of buckle. The minimum value of k s with respect to 
value of wave length is determined from. the. condition 

6k 

"S5T = 0 < C10 > 

which gives (when € does not depend .on wave length) 



(l + 


V 


I o e 

if 3 + =4 - k 0 f 2 -... 


TT 


(Cll) 


Substitution of this value of wave length in equa- 
tion (C9) gives 



+ + S08 ) f 2 - h e2 «i] < C12 > 

The minimum value of k s .with respect to angle of 
buckle is found from the condition 



(C13) 


which gives 


e 2 = 


- /2e , ^ 
f l 3 + 

\ I* 


) - k c f ; 


1/2 


+ fc 


f l(;l + "2 ) - KH 


1/2 


+ 3fg - 


k c f l 


(C14) 
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or> 

60 


If this value of 0 is substituted in (012) , the final 
result is the following interaction formula, in which 
k s is given in terms of k c and the edge restraint e s 


k g 2 = 4C X + 3 Go(C 2 - k c ) + (4C 2 - k c )\j2(2C 1 - Cgk c ) (C15) 


where 



NACA ARR No. L4L14 


21 


REFERENCES 


1. Stowell, Elbridge Z., and Schwarts, Edward 3.: 

Critical Stress for an Infinitely Long Plat Plate 
with Elastically Restrained Edges under Combined 
Shear and Direct Stress. EACA ARR No. 3K13, 1943. 

2. Shanley, P. R., and Ryder, E. I.: Stress Ratios. 

The Answer to the Combined Loading Problem. 
Aviation, vol. 36, no. 6, June 1937, pp . 28, 29, 

43, 66-70. 

3. Anon.': Strength of Aircraft Elements. ANC-5, Army- 

Navy-Civil Committee on Aircraft Design Criteria. 
Revised ed. , Dec. 1942. 

4. Crate, Harold, Batdorf, S. B., and Baab, George VV.: 

The Effect of Internal Pressure on the Buckling 
Stress of Thin-Walled Circular Cylinders under 
Torsion. NACA ARR No. L4E27, 1944. 

5. Hopkins, H. G., and Rao, B. V. S. C.: The Initial 

Buckling of Plat Rectangular Panels under Combined 
Shear and Compression. Rep. No. S.M.E. 3244, 
•British R.A.E., March 1943. 

6.. Wagner, Herbert: Tiber Konstruktions- und 

. Berechnungsf ragen des Blechbaues. Jahrb. .WGL, 

1928, R. Oldenbourg (Miinchen und Berlin), pp. 113- 
125. 

7.. Timoshenko, S.: Theory of Elastic Stability. McGraw- 

Hill Book Co., Inc., 1936, p. 305. 

8. Stowell, Elbridge Z.: Critical Shear Stress of an 

Infinitely Long Plat Plate with Equal Elastic 
Restraints against Rotation along the Parallel 
Edges. NACA ARR No. 3K12, 1943." 

9, Lundquist, Eugene E. , Rossman, Carl A., and Houbolt, 

John C.: A Method for Determining the Column 

Curve from Tests of Columns with Equal Restraints 
against Rotation on the Ends. NACA TN No. 903, 
1943. 



NACA ARR No. L4L14 


22 


TABLE I.- VALUES OP k 0 AND R c 
COMPUTED VALUES 0* k s 


WITH CORRESPONDING 
AND R s 



k s 


Rs 

Exact 

Energy 

Exact 


n 

o 


1.00 

( 0 

<• 2.86 

0 

3.16 

} 1.00 

f 0 

t .535 

.99 

3.09 

3.39 

.99 

.578 

.95 

3.36 

3.67 

.95 

.629 

.90 

. 3.57 

3.88 

.90 

.669 

.80 

3.87 

4.20 

.80 

.725 

* .60 

4.35 

4.66 . 

.60 

.814 

.40 

4.70 

5.04 

.40 

.881 

.50 

4.88 

5.20 

.30 

.913 - 

.20 

5.03 

5.36 

.20 

.942 

.10 

5.19 

5.51 

.10 

.972 

0 

5.34 

5.66 

0 

1.000 

• - 1.00 

6.56 

6.89 

- 1.00 

1.229 

- 2.00 

7.63 

7.91 

- 2.00 

1.428 

- 5.00 

10.22 

10.47 

- 5.00 

1.914 


c = 2 


1.67 

' 1.33 - 

.67 

0 

- 3.34 

- 8.34 

f 0 

l 2.77 

4.07 
5.20 
6.06 
' 9.26 
12.95 

0 

2.96 

4.32 

5.44 

6.30 

9.50 

13.16 

} 1*00 

.80 

.40 

0 

- 2.00 

- 5.00 

f 0 

1 .457 

.672 
.858 
1.000 
1.527 
2.137 


€ = 5 

2.30 

f 0 

l 2.76 

0 

3.03 

j 1.00 

( 0 

' .413 

1.84 

4.27 

4.43 

.80 

.639 

.92 

5.62 

5.81 

.40 

.841 

0 

6.68 

6.88 

0 

1.000 

- 4.59 
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edges elastically restrained against rotation. 







(b) Shear with compression 
3, - Concluded. 
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Figure -4.- Comparison of correct interaction 
curves with a curve formerly proposed for 
infinitely long plates under combined 
shear and compression . 
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Figure 5. - Coordinate system used in Appendix B 
for an infinitely long plate under combined 
shear and transverse compression. 
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Figure 6. - Coordinate system used and wave 
form assumed for energy solution in Appendix C. 
inclined lines indicate nodal lines of buckles. 





